Quantum Friction of Micromechanical Resonators at Low Temperatures 
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Dissipation of micro- and nano-scale mechanical structures is dominated by quantum-mechanical 
tunneling of two-level defects intrinsically present in the system. We find that at high frequencies — 
usually, for smaller, micron-scale structures — a novel mechanism of phonon pumping of two- level 
defects gives rise to weakly temperature-dependent internal friction, Q~ , concomitant to the effects 
observed in recent experiments. Due to their size, comparable to or shorter than the emitted 
phonon wavelength, these structures suffer from superradiance-enhanced dissipation by the collective 
relaxation of a large number of two-level defects contained within the wavelength. 

PACS numbers: 62.65.+k, 62.20.Dc, 62.40.-K 
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Physical properties of physically- and chemically- 
engineered micromechanical systems are of immense fun- 
damental and technical interest . Some of the recent 
spectacular examples of micromechanical structures in- 
clude the measurement of vortex motion in high-T c su- 
perconductors , biomolecular recognition Q , actuation 
of sensors via Casimir force |4| , and shuttling of electron 
charge in a quantum dot 5J. Important to all these ex- 
periments is the oscillation of a particular or a set of 
micron-sized resonators at a resonance frequency deter- 
mined by the geometry and material properties. Changes 
in the "resonant" oscillation frequency or the oscillation 
amplitude mostly determine the magnitude of the force of 
interest, to which the micromechanical structure is cou- 
pled. Detrimental to the detection of force is the damping 
of the resonant structure, quantified by quality factor Q 
or dissipation Q~ l . 

The essential problem inherent to the types of force 
measurement mentioned above is the low quality fac- 
tor Q, observed in small resonators. In many experi- 
ments, intrinsic two-level defects @, Q are found to be 
the dominant cause of internal friction Q _1 in crystalline 
resonators 0, 0, . The observed linear temperature 
dependence of Q~ l oc T || could be explained assum- 
ing reasonable density of defects within a linear response 
theory of two- level systems (TLS) [ill ll^ . Despite the 
success of the TLS mechanism, dissipation at low tem- 
peratures continues to hold many challenges, notably the 
observation of a weak temperature dependence and 
the non-monotonic dependence of Q~ l (T) in silicon res- 
onators 01 • Most recent experiments merely add to the 
list of problems yet to be understood. 

Previous theoretical works on the acoustic response of 
TLS use the linear response approximation or adiabatic 
approximation for nonlinear response [lj, [l^J- The ques- 
tion we would like to explore is whether the anomaly 
can be explained in terms of TLS if both assumptions 
are abandoned. For dissipation in micromechanical res- 
onators, such issues become terribly important, as a new 
parameter space emerges with decreasing size. For ex- 



ample, emitted phonon wavelength becomes much longer 
(3> lpm) compared to the system size at low temper- 
atures, bringing novel correlation effects (lij. extremely 
relavant to mesoscopic structures. 

In this work, we evaluate internal friction Q~ l when 
the two-level defects operate non-linearly and non- 
adiabatically. We show that the nonlinear response of 
Q^ 1 results in a large contribution at low tempera- 
tures at high driving frequencies due to the process of 
phonon pumping. We propose and demonstrate a mech- 
anism of how Q™ 1 can be significantly enhanced by the 
phonon pumping process through cooperative emission 
of phonons. Cooperative emission becomes possible be- 
cause of the long wavelength of the emitted phonons, 
which allows the correlated decay of two-level systems 
contained within the wavelength. Our theory provides a 
mechanism for enhanced dissipation in micromechanical 
resonators within the context of two-level systems. 

The Hamiltonian of a TLS is fully characterized by the 
asymmetry energy A and the tunneling matrix element 
Ao . The asymmetry energy of a TLS is the energy split- 
ting due to the static strain ±^A S and the applied time- 
dependent strain field eo cos Lut: A(t) = A s + 2760 coswi, 
where 7 is the constant of coupling between the acoustic 
wave and the asymmetry energy. Friction due to thermal 
phonons is described by the Hamiltonian for the single 
TLS coupled to phonon fields 0; 
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where V and p are the volume and the mass density of 
the resonator, v a is the sound velocity for polarization a, 
a q,a( a q a ) i s the phonon annihilation(creation) operator. 
7 Q is the coupling constant which relates the change in 
the asymmetry energy to the local elastic strain. 

In crystalline materials with many TLSs, it is known 
that the static asymmetry energy A s has a wide distribu- 
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tion due to the randomness of local strain, while tunnel- 
ing splitting energy A has a well-defined value I 111. For 
simplicity, we consider a rectangular distribution |l2( of 
A s with n number of TLS per unit volume in < A s < 
A™ ax . The underlying assumption, fiuj <C 7e -C A™ aa: , 
holds in most experiments. 

The internal friction Q _1 is related to the energy loss 
per unit volume SE in a cycle of the acoustic wave ; 
Q 1 = ' wnere ~ 2peQ« 2 is the acoustic wave 
energy stored per unit volume |2|. The time-averaged 
phonon emission power P(A S ,T) of a TLS is related to 
Q-h Q-i(T) = ^^-J A ^ P(A s ,T)dA s . Further- 

more, we divide the two-level systems into two kinds 
based on their asymmetry energy A s : (i) A s < 2760, and 
(ii) A s > 2760; which contribute to the total dissipation 
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P(A s ,T)dA s . (3) 



In linear response theory where eo is infinitcsimally small, 
all the TLS belong to category (ii), and Q^ 1 — 0. How- 
ever, even a small experimental value of eo = 10~ 6 will 
give rise to nonlinear behavior 0, which implies that TLS 
with A s < 2-feo may also play an important role[l^|. 

Usually a time-dependent unitary transformation U 
is used for solving time-dependent Hamiltonian in the 
nonlinear regime, so that the TLS has time-dependent 
energy levels given as ±^E(t) JE( t ) = y/ ' A(t) 2 + A 2 ,). 
In the standard adiabatic limit [lj, used to explain 
the low-frequency experiments, the acoustic field fre- 
quency is assumed to be very low so that the additional 
term due to the unitary transformation —ihU^^- — 



'y E(t) 



hui sin cot can be safely neglected. In contrast, 
here we consider the high frequency regime where the 
matrix element of this additional term can not be ne- 
glected in comparison to the diagonal term ^E(t)a z even 
at the level-crossing point E(t) = A : 



> A . 



(4) 



Under this condition, as the tunneling motion between 
the two sites is not fast compared to the rapid time evo- 
lution of the acoustic wave, it is not useful to take the 
TLS bonding states of the two sites as the basis states. 
Instead, it is more natural to consider the localized basis 
states in each local potential minimum and consider the 
tunneling term ^A$<j x as a perturbation. 

We treat the coupling to the acoustic wave adiabati- 
cally and treat the tunneling term and the coupling to 
the thermal phonon pertubatively. The approximation 
is guaranteed by the fact TLS-phonon interaction can be 
treated pertubatively whenever the two-level approxima- 
tion is valid |l8| . Let |Z > and \r > be the left and the 



right localized states respectively. The leading term for 
the transition from one state to the other is the second 
order process involving the tunnelling and emission (ab- 
sorption) of a phonon; at low temperatures single phonon 
process dominates. The thermal averaged relaxation rate 
is written as 
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Here we assume the Debye 



where A = Y\ „ 7? s ■ 
density of states of phonon with a linear dispersion: 
= v a q. 

Time evolution of the occupation number n r (= 1 — 
ni) follows the rate eauation|l4j| for the time-dependent 
asymmetry energy A(i) = 2760 coswt + A s : 



(Jilt ip TJi y Til f 

IT = ~ r(A(t)) ; 



4°\t) = l/(eP*V> + l), (6) 



where j3 — 1/kgT. The typical process is depicted in 
Fig. 1. In the regime A s < 2^€q, the two energy levels 
cross, and at this point population inversion arises. Note 
that the population inversion process ( A — > B — > C ) 
is followed by the spontaneous emission of phonon (C — » 
D). The high-frequency assumption of Eq.Q is crucial 
for the pumping process near the level-crossing point. In 
the conventional adiabatic limit (jeohLO / Ao <C Ao) with 
TLS energies of ±\E[t), level crossing is avoided and the 
TLS in the low (high) energy state remains mostly in the 
low(high) energy state. 

The time- averaged phonon emission power 0] can be 
obtained by solving (EqEJl for n r (t); 
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dt 
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For the general case of asymmetric TLS, numerical in- 
tegration of Eq.lQ gives P as a function of the static 
asymmetry energy A s and temperature T, and Eq.J2J 
subsequently gives Q^ 1 . For the particular case of sym- 
metric TLS (A s = 0), Eq.(JZJ) can be rewritten using an 
analytic solution of n r {t) at T = 0; 
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where r e is a typical relaxation time scale of symmetric 
TLS at zero temperature when wt £ > 1; 



- = 2AA\ 1 t Q . 



(9) 



When A s > 2760, the TLS always stays in the lower en- 
ergy state at zero temperature, hence the TLS does not 
emit a phonon in its steady state, P(A S > 2760, T — 
0) = 0. Since P(A S ,T = 0) is a smoothly-decaying 
function of A s , an analytic formula for Q^ 1 (T = 0) can 
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FIG. 1: The schematic diagram describes the energy-pumping 
process of the two-level atom from the acoustic wave to 
phonons. The crossing curves represent the time evolution 
of the two energy levels with the energy gap of A(i) = 
2760 cos cut + A s . The energy diagrams of the two-level atom 
as a function of the configuration coordinates are depicted in 
the upper panel. 



be obtained using the approximation J^ je ° P(A S ,T = 
0)dA s » 760(^(0,0) + P( 7 eo,0)) = 7 e P(0,0). From 
Eq.JSJl, the zero temperature internal friction Q^ 1 is 
found to be 



Ql\T = 0) 
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(10) 



where Ii(x) is the first kind modified Bessel function of 
order 1. 

Now let us turn to Q^ 1 . For weak strain 2760 < A s , 
where the energy levels do not cross, phonon pumping 
process does not take place. This is the regime where 
the usual linear response theory is a good approxima- 
tion for describing internal friction. In the linear re- 
sponse 

theoryHQjl, Q^CQ in E q-@ for k B T > Tilj > 
H/t (Aq) can be written as[ 



FIG. 2: The internal friction Q~ 
of temperature for u>T t = 1. 
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where E = a/ A 2 + A 2 ,. For A < je < A r s ; 
expression for Q^ 1 is easily found: 
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where f(x) is an almost linear function, defined by 



l/x sinh t ' 

Fig. 2 shows the temperature dependence of Q^ 1 = 
Qi 1 + Q2 1 f° r WT £ = 1) calculated from Eas. l|2l3lt)l7l and 



I12|l . Q^ 1 shows a weak temperature dependence with de- 
creasing temperature, as Qj -1 becomes more important 
than Q^ 1 at low temperatures. The low-temperature 
saturation value of Q^ 1 for ujT e > 1 is well approx- 
imated by ~ Q.3(l/u>T e )(n'y 2 / pv 2 ). Now, let us esti- 
mate r £ from experimental data|8(. Since n^ 2 / pv 2 is 
roughly the high-temperature saturation^ value of 
for the single-crystalline silicon data in Ref.0, we esti- 
mate wy 2 / pv 2 ~ 10~ 6 — 1CP 5 from the experiments at 
uj ~ 10 3 — 10 4 Hz. Since the low temperature satura- 
tion value of Q^ 1 was found to be ~ 10~ 7 — 1CP 6 for 
these frequencies, we estimate wr e ~ 1, which means 
t £ ~ 10~ 3 — 10 _4 sec. However, these values of t € are 
anomalously small compared to the direct estimate us- 
ing Eq.@with reasonable values of parameters^. I?l Hl| . 
By using 7 - leV, e ~ 1(T 6 , A = 10 8 s~ 1 k B 3 into Eq.@, 
we get r e ~ 1 sec. In other words, the low-temperature 
dissipation in experiments is anomalously larger than 
our direct theoretical estimation. 

Here we propose and discuss the origin of the anoma- 
lously small relaxation time T e . We believe that the coop- 
erative emission [l7| of phonons provides a natural mech- 
anism for the short relaxation time. In Dicke superradi- 
ance |l7| , N number of two- level atoms within a distance 
shorter than the emitted photon wavelength emit pho- 
tons cooperatively, therefore the lifetime of the atom in 
the excited state is effectively reduced by a factor of N. 
In the usual treatment of phonon emission of TLS, as we 
have done so far, the emission of phonon is calculated by 
assuming that the TLSs are totally independent. How- 
ever, this assumption is not valid when a number of TLS 
are contained within a distance shorter than the wave- 
length of the emitted phonon. This phonon wavelength, 
\h ~ ^fj (v s is the sound velocity) , can be considerably 
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larger (~ nm), leading to the cooperative emission from 
a large number of two-level systems. 

Now let us consider N two-level systems within a vol- 
ume \p h . In the context of phonon emission, pumping 
of the two-level system is done by the acoustic wave 
applied to the resonator. The many-TLS Hamiltonian 
in Eq.Q can be generalized to H = Hq + H tun + 
H TLS - P h where H = H TLS + H ph = §A(t) J2 l= i °V + 
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The many-TLS state \l, I, I, ...l), where all the TLS are 
in the I state, is the maximally excited state or the ground 
state. We consider the set of many-TLS states generated 
from Htls-pH + H tun , which are 
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where M = -N/2, -N/2 + 1, -N/2 + 2, N/2. These 
states satisfy H TLS (t)\M) = MA(t)\M). As in the case 
of a single TLS, a perturbation analysis finds the transi- 
tion rate wm±i,m from I M) to M ± l): to be 

w M ±i, M (t) = ™(TA(t)) (j =F M\ ± M + (15) 
where 



ib(A) = AA% 



1 - exp(-A/£: B T) 



(16) 



Here, we have used the relation (M ± 1 1 Ei=i a t I -^0 = 
^/(iV/2 =F M)(N/2 ±M + 1). Now the following rate 
equation for the many-TLS remains to be solved numer- 
ically: 



driM 
dt 



WM,M+l(t)n M +l - WM-l,\f{t)riM 

WM,M-i(t)n M -i - w M+liM (t)n M , (17) 



The time-averaged power P^ emitted by N number of 
TLS can be obtained as: 



Pn = ~ — 



2^ 



2-k/u) N/2 

dt MA(t) 

M=-N/2 



dn M (t) 
dt 



(18) 



Note that for phonons independently emitted by the 
TLS, P/v is simply proportional to the number TV; P/v = 
NP. However, in cooperative emission, the lifetime of the 
excited state of a single TLS is N times shorter, hence the 
phonon emission power of a single TLS should be TV times 
larger when uir e > 1 19] . This means that P/v should be 
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FIG. 3: The time-averaged phonon emission power Pjv of N 
symmetric TLSs (A s = 0) at T = 0.1-ye /k B and ujt, = 20, 
which was numerically calculated using Eq.(18). 



proportional to N 2 for wt ( > 1. Our numerical results 
indeed show that the P/v has a quadratic dependence of 
N (see Fig. 3.) 

When N number of TLS are within X p h, one obtains 
Qr 1 by replacing P with its effective value P* = Pn/N. 
The enhanced phonon emission power due to superradi- 
ance P* ~ NP then provides a natural explanation of the 
anomalously- large internal friction at low temperatures. 
Including the phonon superradiance effect, Ea. HlU|) for 
cuTf > 1 can be rewritten as 



- i (T = o) = or 1 ( T = o) 



1 



r«7 



lot* pv*- 



(19) 



where r* « r e /N is the renormalized relaxation time 
and N is the number of cooperating TLSs. Thus, the 
enhancement of Q _1 at low temperature by a factor of 
10 3 — 10 4 compared to our initial estimate in the indepen- 
dent TLS model (Eq.|JTUJ|) indicates that a large number 
of ~ 10 3 — 10 4 TLS emit phonons cooperatively. By tak- 
ing ~ 10/im from A p /j = ^j-, the total number of 
TLSs in the volume A 3 ^ is ~ 10 5 according to Ref . [llT] 
or ~ 10 3 according to Ref.^2], which is reasonably large 
enough number for the cooperative emission by 10 3 TLSs. 

In conclusion, we formulate a theory of internal friction 
of micro- and nano- mechanical resonators, which invokes 
intrinsic two-level defects in the nonlinear regime at high 
frequencies. As temperature decreases, the mechanical 
friction crosses over from the linear regime to the non- 
linear regime, resulting in the saturation behavior as T 
goes to zero. Because of the phonon superradiance, the 
low-temperature saturation value of Q -1 is strongly en- 
hanced by a factor given by the number of two-level sys- 
tems contained within the emitted phonon wavelength. 
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